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We study various solutions of the Sine Gordon model in (1+1) dimensions. We 
use the Hirota method to construct some of them and then show that the wobble, 
discussed in detail by Kalberman [1], is one of such solutions. We concentrate our 
attention on a kink and its bound states with one or two breathers. We study their 
stability and some aspects of their scattering properties on potential wells and on 
fixed boundary conditions. 

1 Introduction 

Topological solitons play an important role in the description of many phenomena in 
physics. In this paper we look at solitons of the simplest model in (1+1) dimensions, 
namely the Sine-Gordon model. 

This model involves a scalar field i^{x, t) and is based on the Lagrangian density given 
by (we set the speed of light to c = 1) 
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This particular model arises in many areas which range from the description of Joseph- 
son junctions [2] to systems with one-dimensional dislocations [5]. The model has also 
been very intensively studied by mathematicians (as it describes spaces with constant 
negative curvature [4j) and by theoretical physicists working in integrable and conformal 
field theories [5]. 

As is well known p] the model possesses kink, antikink and breather solutions. A kink 
solution is a static field configuration which solves the Euler Lagrange equations based 
on ([1]) ie 

ftt - fxx = — —sm{f3Lp) (2) 

and satisfies the boundary conditions Lp{x = — oo) = and = oo) = In. Such a field 
is well known [6] and is given by 

4 

= — ArcTan [exp(m {x — Xq))] (3) 

P 

For the antikink the boundary conditions are interchanged and in the field configuration 
given above there is a '-' sign before {x — xo). 

In addition, the model possesses also so-called 'breather' solutions. These are nonstatic 
field solutions of ([2]) given by 

^/l — uj"^ sin(ma;t) 

io cosh.{m ^/l^^LJ^{x — Xq)) \ 

Here u is a free parameter of the solution which varies from -1 to 1. The breather field, 
which can be thought of as describing a bound state of a kink and an antikink, oscillates 
with frequency mu. 

As the basic Lagrangian is Lorentz covariant - all these field configurations can be 
Lorentz boosted resulting in field configurations moving with velocity v < c = 1. 

Recently there has been some controversy as to whether a kink possesses an internal 
mode [7] - [8]. Such a mode was claimed to exist by Boesch and Willis [9] and then 
disclaimed by Quintero et al [TOj. The mode at stake is a possible internal mode of zero 
frequency which may have arisen in numerical studies of the Sine-Gordon model. Such a 
mode could be a genuine oscillatory mode or a numerical artifact. 

In fact, as is well known [6|, the Sine-Gordon model possesses many solutions in 
addition to above mentioned kinks and breathers. One such solution was recently studied 
extensively by G. Kalberman [1]. He called it 'a wobble' and looked at its properties 
in detail. However, it is not clear from his discussion whose claims his wobble solution 
supports. 

It is worth recalling at this stage that the wobble is only one of many solutions which 
involve a kink with breathers. Hence we have decided to revisit this subject and look in 
some detail at these solutions. In particular, we have decided to present explicit forms 



4 

(y9 = — ArcTan 

P 



2 



of the field configurations which describe some of these solutions, so that they become 
better known outside the 'integrable model community'. We have also decided to look at 
some of their properties, paying particular attention to their stability and their scattering 
properties on defects (here taken in the form of an interaction with potential holes and 
boundaries). This we discuss in the next sections. 



2 The wobble 

The wobble solution of Kalberman [1] involves a field configuration describing a static 
kink and a breather. In his paper Kalberman gives an analytic form of this solution and 
then discusses some of its properties. The analytical form of this solution was obtained 
by Kalberman using the Inverse scattering method of Lamb and Segur. 

In our work we use the Hirota method [TT] of deriving such solutions as discussed in 
detail in the appendix. We derive in section the exact solution describing a kink and a 
breather moving with respect to each other. The solution is given by eq. (H7|) . However, 
the wobble corresponds to a kink and a breather sitting on top of each other and not 
having a relative motion. Therefore, if one sets the the velocities to zero in the general 
solution (H7|) one gets the wobble solution as 

4 \^^^^ sin (ma; t) + ie^"'^ /g-m _^ ^2 ^mVT^x\\ 
= — ArcTan - — . ^ --. — — (5) 



/3 cosh {m \J\ — iO^ + -^^-^^7^ p^emx gjj^ [mujt) 

where w is a frequency varying from —1 to 1, and 
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P = . . n , (6) 



where e = ±1 corresponds to the kink (+1) or anti-kink (—1). 

This agrees with the expression given by Kalberman. As is clear from the field 
configuration depends on one parameter (the frequency of the breather) and so we have 
studied the stability of this field configuration by calculating and its time derivative from 
([5]) and then used the fourth-order Runge Kutta method to simulate the time evolution of 
this configuration (the spatial derivatives were calculated using central differences). Our 
simulations involved looking at a breather-kink system [e = 1) and for j3 and m we took 
(3 = m = 1. The results of our simulations were in complete agreement with the analytical 
expression thus showing that the solution is stable with respect to small perturbations 
(due to the discretisations). 

Next we tried to assess the stability of the wobble with respect to larger perturba- 
tions. We have performed several perturbations, the most important of them being the 
perturbation of the original slope of the kink, {ie in the expression we have replaced 
exp{emx) by exp{Xx) where A 7^ 1). We have performed numerical simulations with 
A = 1.05, 1.15, 1.18, 1.20, 1.24 and 1.3. 
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Figure 1: Energy as a function of time as seen in a simulation started with A = 1.15 

Each perturbation added an extra energy to the system. Such a system was then 
unstable and so it evolved towards a stable wobble emitting some radiation which was 
sent out towards the boundaries of the grid. To prevent the reflections from the boundaries 
we absorbed the energy there. For A close to one - the perturbations were small - hence 
the system returned to its initial configuration (with A = 1.0). For larger values of 
the perturbation the system was more perturbed and often not only kept on sending 
out its excess of energy but also, at regular intervals, altered its frequency of oscillation 
(increasing it) which allowed it to send out even more radiation. In fig 1. we present the 
plots of the time dependence of the total energy as seen in the simulation in which A was 
set at 1.15. 

We have also studied the case of v = Vk = vb 7^ 0. In this case we had a system 
consisting of a kink and a breather moving with a constant velocity. Again, the system 
was stable (we run it with very small v to avoid having problems with the boundaries and 
for larger values of v with fixed boundary conditions) . 

Then we performed a series of simulations in which the initial configuration was sent 
towards a potential hole. This was achieved by making a = ^ in x dependent; i.e. 
we set 

a{x) = 1, for |x| > 5, a{x) = ao < 1 for — 5 < x < 5. (7) 
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Then placing the breather and the kink long way away from the hole {ie from — 5 < x < 5), 
and sending them towards it we could study the effects of their scattering on the hole. We 
have found that the hole can separate the breather from the kink (in one simulation we 
saw the kink being trapped in the hole while the breather bounced off the kink trapped in 
the hole and returned to where the system has originally come from). As studied by two 
of us [12] the scattering of a breather on the hole is very complicated and produces many 
different outcomes; this time we have even more possibilities and so we have decided to 
postpone the further study of this problem to some future work. 

3 Kink and two breathers 

In this section we briefly discuss another interesting solution of the sine-Gordon model; 
namely the solution corresponding to one static kink and two breathers. 

In this shown in the appendix, the field is given, for /? = m = 1, by 

A 

(/? = 4 ArcTan — (8) 



where 



A = -2cotan6'ie^™'^^^ sin(t sin^i)-2cotan6'2e^'°'^^^ sin(t sin^2) + e^ 

- 2 cotan 6*2 4P ^12 ^ e'' cosei+cose^) ^-^ ^-^ ^ ^2 ^2x cose^ 

- 2 cotan Oi (r[p afe ^ ^'^^^ sin (t sin 9i) + pi e^"" 

+ 2a|2ViP2Cotan^icotan6'2e"^e^(=°'^^i+™'^^2) cos (t (sin ^1 + sin ^2)) 

- 2 pi p2 cotan 6*1 cotan 6*2 e' ^ icosei+coso^) ^^-^ _ q^^^ 

+ M+^aypiP2)' e^'e2-'(™^^^+^°^^^) (9) 



and 



Here 



B = 1 + ^ "^"^ + ^ + 2 pi cotan 9i sin {t sin Oi ) 

+ 2p2 cotan 62 sin {t sin 62) 

+ 2ai2^cotan^icotan^2e''('°'^^+^°"^^^ cos (t (sin ^1 + sin ^2)) 

- 2aS^cotan^iCotan6'2e^(=°'^^+™'^^2) cos (t (sin ^1 - sin ^2)) 

I /'^(+) ^(-)A^ „2x (cosei+cosea) 

+ [^^12 0^12 ) e 

+ 2 a[p ^ pi p2 cotan 62 (scosei+cose^) 

+ 2 aiJ) ^ pi pI cotan Oi ""^^^^ sin {t sin Oi) . (10) 

_ 1 - cos Oj (±) _ 1 - cos (6*1 ± 62) , . 

1 + cos tli 1 + cos (6*1 ± U2) 



This solution depends on two constants (6*1, 6*2), which control the frequencies of 
breathers' oscillations. 

We have tested the stability of this solution by using the expression ([8]) to calculate 
(f{t = 0) and ^(t = 0) and then performing a numerical simulation of ([2]). As before, the 
discretisation has produced a small perturbation but the field configuration was stable; ie 
after a long simulation (we run it till t = 5000) the field was indistinguishable from the 
expression and there was no noticeable radiation. Hence we can conclude that this 
field configuration is also stable. 

4 Perturbed Field Configurations 

Given that we have many field configurations which resemble perturbed kinks {ie which 
are given by kinks and breathers) we have tried to see what happens when one perturbs a 
kink and lets it evolve in time. We have looked at various perturbations, paying particular 
attention to configurations which involved adding to a kink an extra perturbation of the 
form 

6^{t = 0) = ,^ , , 6^{t = 0) = — (12) 
cosh(;Ltx) ot cosh(z^xj 

We have looked at various values of A, B, jj, and z/. In all cases the perturbation made the 
kink move and generated many moving breather-like configurations. We let the system 
evolve - absorbing the energy at the boundaries of our grid. This had the effect of slowing 
down the kink. In fig 2 we present the plots of the total energy, and of the potential 
energy, of one such simulation (corresponding to the values A = 0.5, B = 0.5, fj, = 1.0 
and ly = 0.2). The curve in fig 2a shows a steady decrease in total energy down to close 
to the value of the energy of one stationary kink. In fig 3a we show the blowup of the 
total energy at the latter times (at values at which the curve in fig2a may appear to be 
constant). The figures 2b and 3b show the corresponding changes of the potential energy. 

We note some steps of the decrease of the total energy (they correspond to the moments 
when the kink was reflected from the boundaries). The potential has also gradually 
settled as seen from the plot. Its oscillation demonstrates the existence of transient time 
dependent structures [ie breathers). This can be seen from looking at the time dependence 
of individual field configurations. In fig 4. we present the plots of the fields at t=6750 
and t=6753. They show many breather-like structures - the clearest ones being close to 
x=35 and x=-38. 

Thus it is clear to us that a general field configuration will gradually split into moving 
kinks and breathers and some radiation which will quickly move out to the boundaries. 
However, the resultant field configuration is metastable; it still radiates, albeit very slowly, 
and gradually evolves towards a field configuration involving mainly a kink. Whether at 
the end of its evolution we end up with a kink or a kink with some breathers is hard to 
determine. 
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Figure 2: Energies as a function of time as seen when the starting file was given by f|T2l) 
a) total energy, b) potential energy. 
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Figure 3: Blowup of the curves of Fig 2 - to show that even at later times the energies 
are decreasing, a) total energy, b) potential energy. 



7 



-i 



-60 



-40 



-20 



20 



40 



-60 



-40 



-20 



20 



40 



a b 
Figure 4: Field configurations at a) t=6750, b) t=6753. 

5 The energy 

In this section we give a simple formula for the energy of the exact solutions considered 
in this paper. The energy for all the solutions coming from the Hirota method is easy to 
calculate since, as we show below it comes from surface terms. In addition, the energy is 
additive for the non-linear Hirota's superposition of solutions. The Hamiltonian density 
for the sine-Gordon theory ([1]) is given by 



H 



1 1 



sm I -P(p 



-1 2 



(13) 



Replacing the field ip in terms of the Hirota's r-functions as given in the appendix (1271) 

we get 



H = — 



+ 



-)■ 



(14) 



However, this expression can be rewritten as 



H = _(5+-a_)'(lnro + Inn) + 



where 



1 (dlr^ ^ d^TQ ^ din ^ din 



To 



To 



Tl 



Tl 



Tl To 



Tl To 



(15) 



(16) 



8 



The solutions obtained by the Hirota's ansatz (1501) satisfy, besides the Hirota's equations 
(120]) . also the additional equations 



ndiro + Todin - 2 a±ro d^n = (17) 



This fact is proved, in a much more general setting, in ref. [13]. In fact, it is shown 
there that all solutions that can be obtained from the vacuum solution = hj the 
dressing transformations, satisfy (|T7j) . The Hirota's solutions are, of course, of this type. 
Therefore, it turns out that for these solutions one has Hcor = 0, and so 

H = -^a^(lnro + Inn) (18) 

where we have used the fact that (9+ — 9_) = 2 8^ (see (128|1 ). Therefore, the energy 
becomes 

/oo 4 
^dxH = -d,{lnTo + Inn) llZ-oo (19) 

and so is determined entirely by the asymptotic values of n functions. 

In fact, this result, (in a different and less explicit form) has been known for a while for 
kinks and antikinks, to people working in integrable and conformal field theories [IJl [15] . 
Here we have presented it in a form that is more explicit and more easily accessible to 
people working in other areas of physics. For more details of the general proof see [13] . 

Looking at the solutions in the A^-soliton sector, given in (1321) . one observes that the 
asymptotic behavior of the r-function is determined by the exponentials of the FjS given 
in fl3T]) . i.e. 

(20) 



m 

Ti = — 
2 



Zi + — ] X + IZi 



Zi / \ Zr^ 



In f|T9|) we have to evaluate the quantities \^J^^. Therefore, if a given combination 
of exponentials of the FjS dominates the numerator of for x ±oo, then the same 
combination dominates the denominator. Consequently, when the limit x ±oo is 
taken, we are left with the ratio of these two dominant terms which are equal except for 
the constant term (in the numerator) coming from the x-derivative of the exponentials, 
i.e. the terms of the form y (^zi + j-^. Therefore, the contribution to (fT9l) is just the 
terms from the x-derivatives. In addition, tq and ri, given in (1321) . have the same form 
except for the relative minus signs, and so their contributions to (|T9l) are equal, despite 
these minus signs which cancel when the limit x —>■ ±oo is taken. 

Notice that the solution fl32l) contains all possible combinations of the exponentials of 
the FjS, such that each Fj appears at most once. Therefore, it is clear that the dominant 
exponential in the limit x —>■ +oo is exp(X]iFj), where in the sum we include all FjS 
such that Re (^Zi + 7-) > 0. Similarly, the dominant exponential in the limit x — ^ —00 is 

exp (X^i Tj), where the sum involves all F^s such that Re (zi + ^) < 0. Thus, the energy 
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depends only on the modulus of Re (^Zi + j-^ . The parameters Zi can be complex for some 
solutions and so, writing them as Zi = ^-ai+iOi^ gg^g 

Zi -\ = 2 [cos 6i cosh ai — i sin 9i sinh Oj] . (21) 

Thus we conclude that the energy f|T9l) for the solutions fl32|) . is given by 

8 m ^ 

E = —— ^ [ I cos^i I coshttj — i smOi sinho;,] , (22) 

i=i 

and where corresponds to the A^-soliton sector in which the solution lies. The energy 
will then be real for some special choices of the parameters Zi. But whenever this happens, 
the energy is automaticaly positive. Of course, the energy is real whenever the solution 
Lp is real. 

The energies for the solutions we consider in this paper are then given by: 

1. For the 1-soliton constructed in section PV.21 

8m 1 

2. For the breather constructed in section lA.31 



16 m \/T 



2 



-^'breather n-y p, t: ■ (^4) 

p V 1 — 



3. For the wobble constructed in sections [2] and IA.4I 



8m 1 16 m \/l — 

4. For the solution of the kink with two breathers constructed in sections |3] and IA.5I 
(with their velocities set to zero) 

m 16 m r ~ 16 m 



-E'kmk+2breathers — H V ^ ~ ^ "^S^ " ^2 J (26) 



where Ui = sin^j, z = 1, 2. 
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6 Final Remarks 



In this paper we have drawn the attention of the readers to the rich structure of solutions 
(of finite energy) of the Sine-Gordon model. Even in the one kink sector there are solutions 
involving, in addition to the kink, also many breathers. As the energy of each breather 
depends on its frequency (and vanishes in the limit of this frequency going to 1) the 
extra energy, due to these extra breathers, does not have to be very large. The solutions 
appear to be stable and this stability is guaranteed by the integrability of the model. 
We have tested this numerically and have found that small perturbations, due to the 
discretisations, do not alter this stability. To change it we need something more drastic 
- like the absorption or the space variation of the potential {ie the coefficient of the sin 
term in the Lagrangian). But even then the effects are not very large - one sees splitting 
of breathers etc but no 'global annihilation'. 

At the same time we have looked at the total energy and the energy density of a 
general solution. We point out that the total energy is determined by the asymptotic 
values of the fields {ie is given by (1191) ). This result has been known earlier (for kinks 
and antikinks) to people working in integrable and conformal field theories but may not 
be known generally. We have given a general proof of this in a separate, rather technical, 
paper [13] but we have also checked it explicitly for all field configurations involving up 
to 5 kinks (or antikinks). 

Our results do not answer, definitively, the question as to whether a kink possesses 
an internal mode of oscillation or not. The perturbation of the kink we performed in 
section H] did not produce any oscillatory internal mode. On the contrary, all the energy 
given to the kink by the perturbation was used to produce breather like excitations which 
died away very slowly. In fact, the extremely slow decay of these excitations indicate the 
difficulty of settling down the issue of the existence of the internal mode. If simulations 
are not done very carefully and runned for very long times, then this fact can lead to 
incorrect interpretations. 

As we have pointed out, one can get oscillitory kink configurations by constructing 
exact solutions corresponding to the stationary supperposition of a kink and one or many 
breathers. Such a case of a kink and a breather was named 'a wobble' by Kalberman [1]. 
In the case of the wobble the frequency of oscillation cannot be greater than 1, and the 
energy of the oscillation goes to zero as the frequency approaches 1. Boesch and Willis [9J 
claimed to have seen an oscillatory mode of the kink just above the phonon band, i.e. just 
above 1. If that is so, the wobble does not correspond to that mode. It is true however, 
that the frequency of the wobble can go above 1 by a Lorentz boost. One has then to 
settle the issue of how precise the simulations of Boesch and Willis were to separate this 
effect. If one considers the exact stationary supperposition of a kink and two breathers 
one can get frequencies of oscillations greater than 1, as shown in ([8]). However, this is in 
a context different from that discussed in the literature where the considered frequencies 
are studied in the linear approximation. In addition, although our simulations have shown 
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that the kink plus breathers are stable solutions against the discretisation, they can be 
pulled apart by scattering through a hole, as mentioned at the end of section [21 The 
fact that a kink solution of the equations of motion can involve many breathers makes 
the problem of the zero mode difficult to resolve analytically (and in the literature it is 
discussed only in the 'linear approximation') and almost impossible to resolve numerically. 
The rich spectrum of the solutions and the appearance of many breathers makes this task 
particularly hard to perform. It would be interesting to see whether these extra breathers 
play a significant role in any physical applications of the model. 

A Appendix 

A.l Hirota's method 

In order to solve the sine-Gordon equation ([2]) by the Hirota's method we introduce the 
Hirota's r-f unction as 

^ = ^log^. (27) 

fj To 

Using light-cone coordinates 



x± = -{t±x) d± = dt± 

and replacing fl27j) into ([2]) we get 



= -dl = d+d^ 



:28) 



To 



t2 



1 



Ti rf 

Imposing both sides to vanish we get the Hirota's equations for the sine-Gordon model: 

,2 



To d+d^TQ - a+To a_To + — (tq - tI 

The solutions of (12^ are obtained by the Hirota's ansatz 



0, 



0. 



(29) 



N 



= 1 + £ E t'' e^' + e'Y. + E 



N 



{")^ri+r,- 



N 



i,{")pr,+r,+rfc 



+ 



1=1 



i,i,k=l 



where 



m \ ZiX+ — 



a = 0,1, (30) 



(31) 
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and where Zi are arbitrary (complex) parameters. Replacing fl30p into fl29l) and expanding 
in powers of e, one obtains the coefficients 6*^"'''s recursively. The series in (I3UI) truncates 
at order leading to an exact solution, given by 

N N / \ ^ 

1=1 h<h=i V^'i ^'2 / 

Z.<.^/3 = lV^'i+^j U^i+^j 

••• + (-ir^ n ( '"''''A llaie^^^^^ for « = 0, 1, (32) 

fci<fe2 = l V^fcl + ^'=2/ i=l 

where Oj are arbitrary parameters. 

A. 2 The one-soliton solutions 

The one-soliton solutions correspond to = 1 and the following choice of parameters in 
(I32D 

ai=ze-=^'^° zi=ee-'^ e = ±1 a real. (33) 

Then the argument of the exponentials become 

TTl 

Ti = e'y{x — vt) 7 = m cosh a = t) = tanha. (34) 



Therefore, 1^ and ([32]) gives 

4 

= — ArcTan exp [e 7 (x — v t — Xq)] . (35) 
A. 3 The breather solution 

The breather solution lies in the 2-soliton sector, and it is obtained by taking N = 2, and 
the following choice of parameters in fl32|) 

z, = 6-"+^'^ Z2 = zl ai = i a2 = -al. (36) 

tanfc' 

Then we have that r2 = and introducing Fi = Fi + 77 + we obtain 

fi = Ffi + «F, (37) 

with 

^ m cos 6 , , ^ m sin 9 , . ^ , . 

Tr = ^== ix-vt)+r] Ti= (t-vx)+^ (38) 
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and 
Thus 

If one now takes 



V = tanh a. 



4 (cotan6') cosF/ 

= - ArcTan — . 

p cosh i ji 



(39) 
(40) 



one gets 



1] = V = ^ 



= — ArcTan 



TV 



UJ 



= sin 9 



- <e < - 

2 - - 2 



sin {mut) 



UJ 



cosh {m Vl — 



(41) 



A. 4 The wobble or the kink with a breather 



To have a configuration describing a kink with a breather we take the following values for 
the parameters of the 3-soliton solution (A^ = 3) in ( l32l) 



ai 



-ag+tti 



Z2 = Zl 

a2 = —al 



Z3 = ee 



-an 



(42) 



tan 6* 
Then one gets 

^3 = £ JK — vk t) + tjk Jk = cosh ax vk = tanh ax (43) 

and again one has r2 = F^. Then introducing Fi = Fi + r^^ + i^s, one obtains 

fi = F^ + zF, (44) 

with 

Tr = 



m 



m 



COS 9 {x — Vst) + TjB 

sin 6 {t — vb x) + 



and vb = tanh as . Next we define 



Zl - Z3 
Z\ + Z3 



pe 



i <f> 



with 



cosh (aB — dK) — s cosO 

P = — n ^ n ^^'^ 

cosh [aB — oik) + £ cos 9 



2 ArcTan 



e sin 9 



sinh {aB — ax) 



(45) 



(46) 
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Then 



If one now takes 
and denotes 
then 



4 [2 (cotan6') cosT/ + e"^^ ie'^^ + e^A 

P [(e-rfl + 6^^) - 2 (cotan6') pe^^ cos [Ti + 



(47) 



Vk = Vb = vb = vk = 



= sm9 <e < - 

2 - - 2 



IT 



P = 



and so 



_ 1 -eVT^^ 
1 + 



= — ArcTan - — j \— 

P cosh (m a/I — o;^ xj + ^^"^ pe^'"^ sin (mcu 



t) 



(48) 



(49) 



A. 5 The kink with two breathers 

To get a field configuration describing a kink and two breathers that are all at rest and 
located at the same position we take the parameters of the 5-soliton solution (A^ = 5) in 
^ as 



zi = e* 



take the parameters 

402 -y — -y* 
Z4 Zg 



Z2 = Z3 = e 



Z5 = £ 



and 



Oi = —02 = — cotan 9i 



Then setting m = (3 = 1, one gets 
Fi = = X cos 61 + it sin 9i 
In addition, we have 



as = —04 = —cotan ^2 



as = i. 



(50) 
(51) 



F3 = F4 = X cos ^2 + it sin6'2 F5 = e x. (52) 




tan^^i, 

VZ4 + Z4 
' Z2 — -25 \ ^ 1 — £ COS 6*1 _ 



.2:2 + 

^ /Z4j-Z5\2 ^ ^ 

V-24 + Z5J 
_ f Z2 - Zi \^ _ _ 

V Z2 + -24 / 



2^3 — ^4 \ ^ , 2/1 

- — — ) = - tan^ 62, 
= Pi, 



1 + e cos 6*1 

1 — £ COS 6*2 _ 

1 + e cos fc'2 

1 - cos {6, - ^2) _ (-) 

1 + cos (01-^^2) ' 

1 - COS (01 + 6*2) _ (+) 
— = <^12 • 



-21 +-24/ V-22 + -23/ 1 + COS (6^1 + 6'; 



(53) 
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The r-functions now become {a = 0,1] 



,e X 



To, = I 

+ {-IT I [-2cotaneie^ sin (t sin ^i) - 2 cotan^a sin (t sin02)+e 

^ e^^™'^^ +e2^^°'^2 + 2picotan^ie"'e^^°'''i sin (t sin^i) 

+ 2 p2 cotan 6*2 ^ sin {t sin 02 ) 

+ 2 afe^ cotan ^1 cotan 6*2 6^^™'^^^+™'^''^) cos (t (sin^i + sin6'2)) 

- 2 a[p cotan cotan 02 cos (t (sin Oi - sin ^2)) 

+ {-iTi [-2cotan6'2 aj2^ai2^e^(2^°'^''i+'=°'^''2) sin (t sin 6'2) + e^^ ""^^^ e"^ 

- 2 cotan ^1 aj^^ ^{2^ (=°'^^i+2 ^"'^e^) gin {t sin6'i) + p^ e^^ ^"'^^^ e"^ 

+ 2 a{2ViP2 cotan 6*1 cotan 6*2 e"^e^(™'^''i+™'^''^) cos (t (sin 6*1 + sin ^2)) 

- 2 pi p2 cotan Oi cotan ^2 ^ ^1+^°' cos (t (sin - sin ^2) ) 

I „2x (cos ei+cos6»2) 

+ 2a;j^ai2V?P2Cotan^2e""e"(2-°^^^+=-^2) sin (t sin ^2) 

+ 2ai^^ai2ViP2Cotan^ie""e"(^°^^i+2"°^^^) sin (t sin^i) 

+ (-1)^2 {(^Y2^ (^vi' Pi 92^ e"^e2^(™'^^+™'^2). (54) 
This gives us the expression mentioned in section [31 
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